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Abstract 

We calculate the electric quadrupole to magnetic dipole transition ratio 
E2/M1 for the reaction 7^ — ► A(1232) in the chiral quark soliton model. 
The calculated E2/M1 ratio is in a good agreement with the very new exper- 
imental data. We obtain non-zero negative value for the electric quadrupole 
N — A transition moment, which suggests an oblate deformed charge struc- 
ture of the nucleon or/and the delta isobar. Other observables related to this 
quantity, namely the N — A mass splitting, the isovector charge radius, and 
isovector magnetic moment, are properly reproduced as well. 
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The ratio of electric quadrupole to magnetic dipole amplitude (E2/M1) for the reaction 
7 + N — > A(1232) is a quantity sensitive to a presence of charge deformations in the baryon 
structure. The most reliable phenomenological estimate E2/M1 = (—1.5 ± 0.2)% so far 
comes from detailed analysis [|lj of the available photoproduction data, assuming the most 
general 7iVA gauge coupling and taking into account the unitary condition via Watson the- 
orem. Very recently a very precise 7r°( + )-photoproduction experiment || has been performed 
at MAMI, Mainz, which allows for a direct model-independent estimate of the ratio E2/M1. 
The preliminary result[| E2/M1 = (—2.4 ± 0.2)% || confirms the negative sign suggested 
from the analysis |T| and shows a larger £'2/Ml-asymmetry. This non-zero negative value 
is a clear indication for the presence of an oblate charge deformation in the nucleon or/and 
delta and as such it imposes strong constraints for the effective models of baryon struc- 
ture. Wirzba and Weise have investigated the E2/M1 ratio in a modified Skyrme model 
which includes stabilizing fourth- and sixth-order terms ||. They obtained values between 
—5% and —2.5% depending on the coupling parameters of stabilizing terms. However, the 
other related observables, namely the N — A mass difference, charge radii, and the isovector 
magnetic moment, are not properly described. 

In the present work we study E2/M1 ratio for the process 7 + iV —>■ A (1232) in the chiral 
quark soliton model (for review see ref. [|]]). We employ the simplest SU (2)- version of the 
model with up and down quark degenerated in mass, which is based on the semibosonized 
Nambu Jona-Lasinio lagrangean HU : 

£ = ^(-i 7 ^+m + M[/ 75 )*, (1) 

with auxiliary meson fields 

U(x) = e^'W* , (2) 

constrained on the chiral circle. The model is non-renormalizable and a finite cutoff is 
needed. The latter is treated as a parameter of the model and together with the current 
quark mass mo is fixed in the mesonic sector to reproduce the physical pion mass m n and 
the pion decay constant f n . The last model parameter, the constituent quark mass M, can 
be related to the empirical value of the quark condensate but actually it still leaves a broad 
range for M.| 

In the model the baryons appear as a bound state of N c (number of colors) valence quarks 
coupled to the polarized Dirac sea. Since the model lacks confinement the proper way to 
describe the nucleon is to consider J7[ a correlation function of two iV c -quark currents with 
nucleon quantum numbers JJ^,TT^ at large euclidean time-separation: 

lim U N (T) = (J N (x,+T/2)4(y,-T/2)) 

1 — ► OO 

1 r Nc 1 

^V { llrrJ { lZn, J VUj[(T/2,S\^\ - T/2, y) fi9i e N ^ D ^ , (3) 



In fact, the authors have faced this experimental result after completing the calculations. 

2 Actually, in order to obtain a good description of the baryonic properties the mass M has to be 
chosen around 420 MeV (see review [||). 
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where current Jn is a composite quark operator ||: 

MZt) = j^j rS.TTa *M(^t)-*w^0- ( 4 ) 

In the above path integral the quark fields are integrated out. The effective action includes 
the Dirac operator 

D(U)=idt-h(U), (5) 
with one-particle hamiltonian given by 

h(U) = + /3MTP 6 + (3m . (6) 
i 

In the model the nucleonic solution is obtained in two steps. In the first step, in leading 
order in N c the integral over the meson fields U in eq.((^) is done in a saddle point approx- 
imation. To that end we look for a stationary localized meson configuration (soliton) of 
hedgehog structure 

U(x) = e l ™ p(x \ (7) 

which minimizes the effective action. Actually, the soliton solution is found by solving the 
corresponding equations of motion in an iterative self-consistent procedure 0. Since the 
hedgehog soliton U(x) does not preserve the spin and isospin, as a next step we make use 
of the rotational zero modes to quantize it. It is done assuming a rotating meson hedgehog 
fields of the form 

U(x,t) = R(t)U(x)R + (t) (8) 

with R(t) being a time-dependent rotation SU(2) matrix in the isospin space. It is easy to 
see that for such an ansatz one can transform the effective action 

Tr log £>([/) = Trlog(D(C7) - Q) (9) 

in order to separate the angular velocity matrix: 

n = -iR + (t)R(t) = ^n a r a . (io) 

Since the angular velocity is quantized according to the canonical quantization rule, it ap- 
pears as Q a ~ This allows one to consider Q as perturbation and to evaluate any 
observable as a perturbation series in fl which is actually an expansion in The path in- 
tegral over R(t), which appears due to the ansatz (gj), determines the spin-flavor structure of 
the nucleonic solution. For given spin J, J 3 and isospin T, T 3 this structure can be expressed 
through the Wigner D function 

\N,T 3 J 3 )(R) = (-1) t+t *V2TTTdT=J J3 (R). (11) 

Very recently, the nucleon electromagnetic form factors have been calculated || in this 
semiclassical quantization scheme up to the next to leading order in angular velocity. The 



3 



results for the constituent quark mass M ~ 420 MeV are in fairly good agreement with the 
experimental data. It should be also noted that the 1/N C rotational corrections improve 
considerably the theoretical value for the isovector magnetic moment [110] . 

Obviously it is worthily to consider in this model also the ratio E2/M1 for the process 
7 N — ► A whose value is well settled by the recent experiment in Mainz |2| . It is defined |TT 
as a ratio 



E2 _ 1 M E2 (k, A = +1 ; p( J 3 = -|) -> A+( J 3 = + f )) 
Ml " 3 M ^(k, A = +1 ; p(J 3 = -1) - A+(J 3 = +§)) 



(12) 



of electric quadrupole amplitude M £2 to magnetic dipole one M Ml . Here is the momentum 
of a photon of helicity A in the A rest frame: 

_ Mj-M% 

k ~ 2M A ■ {l6) 
Both amplitudes are related to the corresponding matrix element of the isovector current: 

j$(x) = ^(x)Yj^(x). (14) 

For the electric quadrupole amplitude M E2 , according to the Siegert's theorem (see for 
instance [I2|), one can use the zero component j® as well the space component - V • j a . 



Similar to ref. || we decide to express the amplitude M m via the N — A transition matrix 
element of j® (charge density) of the current j M : 

M E2 (k, X;N^A) = Vl5^ J d 3 x (A\jl(x)\N)Y 2X (x)j 2 (kx), (15) 

The reason is that in the present model this quantity can be calculated || directly in terms 
of quark matrix elements whereas the for V • j a one should use supplementary the classical 
equations of motion (saddle point), for which the ansatz (|j) is apparently not a solution. 
The amplitude M Ml is directly related to the N — A transition matrix element of the 



space components j*: 



M M \k, A; N - A) = -A^ J d 3 x(A\(x x j 3 ) x \N) 3l (kx). (16) 
In the case of A = +1, Jf = — | = +|, we have 

M £2 = +^ / drr 2 j 2 (kr)p E N l(r) , (17) 

and 

M Ml = -3 J drr 2 Jl (fcr)p^i(r) , (18) 

respectively. 

For the matrix element of the isovector current j£ we follow the line of ref. || and in fact, 
we use the results presented there. Here we will only sketch the derivation. We evaluate the 
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expectation value of the quark bilinear operators, \Ed7°7 M Tp-\I/, represented by the euclidean 
functional integral [|7| with lagrangean (JJ): 

(N'J\^(0)^^(0)\N,p) = lim \ [ d 3 xd 3 yS T '- p4T - if3 ' +i ^ 

2 T'— >+oo £ J 

x Jvujw J V^J N ,(x',V) ^( ) 7 V^(0)4(^T)e-/ d4 ^ tD(c/) * . (19) 

Integrating out the quarks in (|19|) it is easy to see that the result is naturally split in valence 
and sea parts. After that we follow the same steps as for correlation function (Q). First, 
we integrate over the meson fields in saddle point approximation (large N c limit). As a 
second step, we take into account the rotational zero modes using the ansatz (||). Due to 
the collective path integral over R(t) our scheme^] involves a time-ordered product of the 
collective operators 

n a (R(t)) = -iTr(R?(t)R(t)r a ) and D ab (R(t)) = ^Tr(R^t)r a R(t)r b ) , (20) 

which appear after the expansion in Q. 

In the above scheme the matrix element of the time-component of the isovector current, 
calculated in the semiclassical quantization scheme, includes only terms linear in Q ~ 1/A C 

(^\m\N) = ^{^nUe m ,e n )<^l(x)r a ^)) (n\r e \m) 

E ^(^)T $ W (f))M|T>)}(A,J^|{J c ,L>3a}|iV,J3T3), (21) 

whereas the matrix element of the space-components of the isovector current includes leading 
order terms ~ Q° as well as next to leading order ones ~ Q (1/N C ): 

(A \j k 3 (x)\ N) = i\J ($L(x) 7°7*r a ® val (x)) (A, J' 3 T> \D 3b \ N, J 3 T 3 ) 

i f$L ; (x)7°7 fc r b $ n (x) N ) (n |r c | vol) 

+ E - \ A,J>T>\{J c ,D 3b ]\N,J 3 n) 

-E<(^)(^(^)7Vr^ n (f)) (A, J' 3 Ti \D 3b \ A, J 3 T 3 ) 

n 

+5 E n M2^ e n )($| n (x)7 7V & <I>„(x)) (n \r c \ m)(A, |[J C , D 3b }\ A, J 3 T 3 )| . (22) 



Here is the moment of inertia, and $ n and e n are the eigenfunctions and the eigenvalues 
of the hamiltonian (H). The regularization functions TZq, 7Zmi an< ^ ^M2 can be found in ref. 



The details of this procedure can be found in ref. 
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For completeness we present the final results for the electric quadrupole density pfjj\{r), 
split in valence and Dirac sea parts, 

PnT 1 (t) = £ (^(r)||[r 2 ®r«]«||$ n (r))(m/||r«||n), (23) 



and 



P^>) = f^f E ^(^,e m )($ n (r)||[F 2 ®r( 1 )]( 1 )||$ m (r))(n||r( 1 )||m) (24) 

n,m=all 



as well as for the magnetic density Pna(x): 



^(r) = -A^(<Mr)||75p (1) ® ® rW](°)||^(r)J , (25) 



P^V) 



E <i(e„) V / 2^n + l($n(r)||7 5 [[r {1) a]« ® ||* B (r) 



n=all 



(26) 



p M r ^) (r) = _ 



^4 E ^^(^MiN^^ 
6 36 n+val e n - e val V J 

(27) 



p Mi ; ^ )(r) = _^ i E <2(£n _ em) ^ n(r) || 75[[f (D x ,(!)](!) x r W]W||$ m (r)) 

n,m=all 

x(n||r (1) ||m) . (28) 

On the other hand, using (0) and dl6|) in the approximation k ■ R <^ 1, where i? is the 
nucleon charge radius, we get following simple formulae 

M E2 = -^=fc 2 (Q„) Ka , (29) 

MM1 = <30) 

where (Q zz )na is the electric quadrupole transition moment and /xjva is the transition mag- 
netic moment: 

Una = "~y|^ / = 1 ■ (^1) 

Here \ii=\ = p p — \i n is the isovector magnetic moment. Using fl29|) and fl3"0|), in the k- i?< 1 
approximation one can relate the ratio E2/M1 to the electric quadrupole iVA transition 
moment: 
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— = -kM N ^^. (32) 
Ml 2 /i M K J 

It should be noted that in contrast to the Skyrme model, (Q ZZ ) NA cannot be directly related 
to the isovector charge radius. The corresponding transition charge density Pna( x ) has a 
more complicated structure, including a spherical harmonics tensor Y^, which acts on the 
quark wave function as a projector for the charge deformation. 



In the numerical computations we use the method of Ripka and Kahana [13| for solving 
the eigenvalue problem in a finite quasi-discrete basis. 

In table |, we present our results for the ratio E2/M1 as well as for some related observ- 
ables, namely the isovector charge m.s. radius, the N — A transition magnetic moment Una, 
the N — A mass difference, and the quadrupole electric transition moment (Q zz )na, for three 
different values of the constituent quark mass M. We compare our results with the exper- 
iment as well as with the numbers of ref. ||. With constituent quark mass M around 420 
MeV we obtain for the E2/M1 ratio values between —2.5% and —2.3% quite in agreement 
with the last experiment data . It should be mentioned that for the same values of the 
constituent quark mass M the nucleon properties (including also the nucleon form factors) 
are reproduced fairly well ||. The only exception is the N — A transition magnetic moment 
which is underestimated by 25%. Our results for other observables in table | show an overall 
good agreement with the experiment which, however, is not the case for the Skyrme model 
calculations ||. The Skyrme model results show a strong underestimation of the isovector 
charge radius and N — A mass splitting whereas the isovector magnetic moment is strongly 
overestimated. 

In the table [I] we also present the results for the ratio E2/M1 in the k ■ R <C 1 ap- 
proximation. Despite that this approximation is not justified it seems that it works in 
practice satisfactorily: the numbers using the formula fl3"2] ) overestimate the exact results 
by not more than 20%. From the relation (|32|) we also get a rough estimate for the elec- 
tric quadrupole transition moment {Q zz )na = —0.026 using the experimental values for 
E2/M1 = —2.4 ± 0.2 |2| and ^na — 3.3, which is not far from our model prediction —0.02. 
This negative non-zero value indicates a presence of an oblate type of charge deformations in 
the nucleon or/and delta structure. It is interesting to mention that the dominant contribu- 
tion to the (Qzz) na in the NJL model comes from the Dirac sea. It can be seen in the table | 
as well on the figure [JL where the electric quadrupole transition moment density, separated 
in valence and Dirac sea parts, is shown. It means that the main charge deformation is due 
to the polarized Dirac sea, whereas the valence quarks are almost spherically distributed. 
Since using the gradient expansion, the polarization of the Dirac sea can be expressed in 
terms of the dynamical pion field - pion cloud, one can think of the nucleon or/and of the 
delta as consisting of an almost spherical valence quark core surrounded by a deformed pion 
cloud. 

In summary, we study the electric quadrupole to magnetic dipole transition ratio E2/M1 
for the reaction 7iV — > A(1232) in the chiral quark soliton model. The calculated E2/M1 
ratio for the constituent mass around 420 MeV is in a good agreement with the very new 
experimental data. We obtain a non-zero negative value for the electric quadrupole transition 
moment, which suggests an oblate deformed charge structure of the nucleon or/and for the 
delta isobar. Other related observables, namely the N — A mass difference, the isovector 
charge radius, and the N — A transition magnetic moment, are properly reproduced as well. 
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FIGURES 

FIG. 1. Electric quadrupole N — A transition moment densitiy, separated in valence and Dirac 
sea parts. 
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TABLES 



TABLE I. Ratio E2/M1 and some related observables, calculated in the NJL model for three 
different values of the constituent massM = 400, 420 and 450 MeV, compared with experimental 
values. The Skyrme model results p| are also presented. 



Constituent quark mass M 



Quantity 


400 MeV 


420 MeV 


450 MeV 


Skyrme 


Exper. 


total 


sea 


total 


sea 


total 


sea 


<r 2 > /=1 [fm 2 ] 


0.88 


0.35 


0.84 


0.37 


0.79 


0.41 


0.60 


0.86 


HAN [n.m] 


2.34 


0.57 


2.28 


0.58 


2.20 


0.60 


6.18 


3.33 


M A - M N [MeV] 


255 




278 




311 




199 


294 


(Qzz)an [fm 2 ] 


-0.020 


-0.014 


-0.020 


-0.015 


-0.021 


-0.016 


-0.028 


-0.026Q 


M m 


0.012 


0.008 


0.013 


0.008 


0.013 


0.009 






M M1 


-0.189 


-0.041 


-0.186 


-0.042 


-0.182 


-0.043 






E2/MI [%} 


-2.19 




-2.28 




-2.42 




-2.6 


-2.4±0.2 


E2/M1 [%] 4 


-2.64 




-2.79 




-2.99 






-2.4±0.2 



4 Using eq.(H) 
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1 1 1 1 1 1 1 1 1 1 r 

1 2 3 4 5 

Radius [fm] 



1 



